Естественные науки 


УДК 536.46 

НАХОЖДЕНИЕ УСЛОВИЙ ДИФФУЗИОННО-ТЕПЛОВОЙ УСТОЙЧИВОСТИ ПЛАМЕНИ 
С ПРИМЕНЕНИЕМ МОДЕЛЬНОЙ ФУНКЦИИ СКОРОСТИ ХИМИЧЕСКОЙ РЕАКЦИИ 

К. О. Сабденов 

Томский политехнический университет. 

Е-таіІ: $аМепоѵ@к21.р!іШри.ги 

Решается задача о диффузионно-тепловой неустойчивости ламинарного пламени с ненулевой толщиной зоны химической ре- 
акции на основе модельной функции источника. Предположение переменности скорости пламени приводит к формуле 
Маркштейна, связывающей скорость движения фронта горения с кривизной этого фронта и практически к абсолютной неустой- 
чивости горения. При постоянной же скорости пламени и числе Льюиса Ее > 1, реализуется только апериодическая потеря устой- 
чивости, а при Ее < 7 — только периодическая. Если энергия активации в единицах произведения газовой постоянной на темпе- 
ратуру пламени выше 6, то абсолютная устойчивость по отношению к возмущениям с любыми длинами волн возможна только 
в малой окрестности Ее = 1 и Ее = 0. 


В наиболее общей постановке задачи на иссле- 
дование диффузионно-тепловой неустойчивости 
рассматривается поведение горения по отношению 
к произвольным возмущениям, приводящим к 
пространственному искривлению плоского лами- 
нарного пламени, движущегося относительно ис- 
ходной горючей смеси со скоростью ѵ„, которую 
называют нормальной скоростью распростране- 
ния. Исследованию неустойчивости ламинарного 
пламени посвящено большое число работ [1-6 и 
др.]. Первое математическое исследование т.н. од- 
номерной устойчивости проведено Д. Розеном [1], 
Г.И. Баренблатгом, Я. Б. Зельдовичем, Я. И. Кане- 
лем [2, 3]. Более поздние [5, 7 и др.] исследования 
показали, что область устойчивого горения опреде- 
ляется не только величиной числа Льюиса, но и па- 
раметром ф=Е(Т ь -Т 0 )/2КТь»1, характеризующим 
скорость химической реакции. При достаточно 
больших значениях данного параметра область ус- 
тойчивого горения локализуется в малой окрест- 
ности вблизи прямой Бе = 1 на плоскости (Бе, у/), 
являющейся единственным условием устойчивос- 
ти в пределе у/-»со [7, 8]. Во всех вышеуказанных 
работах скорость химической реакции моделирует- 
ся 5-функцией, что приводит к несовпадениям ре- 
зультатов различных авторов [4, 7]. Поэтому предс- 
тавляет интерес исследование с модельными функ- 
циями скорости реакции, более близкими к закону 
Аррениуса. 

Б Плоский фронт пламени. В простом случае 
брутто-реакции первого порядка скорость Жпроте- 
кания химической реакции описывается формулой 

Ж = Ш 0 ех р(--^). (П 

Аналитическое исследование задач нестационар- 
ного горения с применением (1) является очень 
трудоемкой работой. Поэтому вместо этого выра- 
жения примем более упрощенную форму 

Ж=М ое х Ѵ \--^уТ-Т,) ( 2 ) 

с температурой 7 . воспламенения, вид которой оп- 
ределяется из требования совпадения скорости 


плоского пламени с аналогичным выражением из 
теории Зельдовича-Франк-Каменецкого при боль- 
ших энергиях активации. Скорость реакции (2) 
сохраняет важнейшие свойства закона Аррениуса - 
его сильную зависимость от температуры и нели- 
нейность. Формула (2) была успешно применена в 
работе [9] при аналитическом расчете скорости 
движения плоского стационарного фронта пламе- 
ни при произвольных числах Льюиса. 

С использованием (2) плоский стационарный 
фронт пламени в газовой смеси описывается систе- 
мой уравнений 


ит-т.), 


кт 


ш . 


Ы N 


ѵ = 

" ах’ сіх 11 


-Ж 0 ехр| --^г|77 (Т-Т). 


( 3 ) 


С введением безразмерных параметров и масш- 
таба ѵ. скорости 
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система (3) принимает вид 

м>°^ = ^г + ]У, \ѵ°ф-=Т е^-+Ж, 
ах СІХ Фх ах 

Ж = а(\-Ъ)г](и-и,\ а = Щ г ехр[--^- |, 

V" V Лі, 


(4) 


где іѵ 0 - безразмерная скорость плоского стацио- 
нарного пламени. В дальнейшем верхний индекс О 
будет означать принадлежность символа к такому 
пламени. 

Задаче о распространении пламени по уравне- 
ниям (4) соответствуют граничные условия 

х — > -со : и 0 = Ь° = 0; х — > +со : Фи° /Ох = аЬ° Ых = 0. (5) 


Помещая место разрыва (2) в точку х = 0 и прис- 
ваивая индексы 1 и 2 температуре и° и выгоранию Ь° 
соответственно при х < 0 и х > О, приведем распре- 
деления м°(х), Ь°(х): 
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х < 0 : и° = и , ехр(ѵѵ°х) = 


Ь" ==■] 1 — ЖА. | е хр(ѵЛс/Ье) = 


к + \ѵ° 
А- 2 Ье 


х > 0 : и° = 1 — 


ехр(тѵ°х), 

ехр(тѵ°х/Ье), 
\ѵ°к < 


к + \ѵ° 
л/(тѵ 0 ) 2 +4яЬе - \ѵ 


ехр(-Ах), Ь 2 ~ 1 — ~ — ~ехр( Ах), 


к п к 


2Ье 


-ехр - 


7? Т, 


(тѵ 0 ) 2 = 


1-м, 


Ье + (1-м,)/и, ’ 


д$/дІ = Щ, А'$ = -Х 2 р, Я 2 = А 2 +А, 2 . 

Параметр в стационарном режиме горения яв- 
ляется постоянной Маркштейна. 

Подстановка и, Ь, ѵѵ в (7) дальнейшая линеари- 
зация и использование (2) в промежуточных вык- 
ладках после несложных преобразований приводят 
к уравнениям для Р, О 


(6) аСу_ар_ {П + х 2^ р + Ш_Р = _Ш' 0+? іфІ_ 

сіх сЬс ди° дЪ" Лх 


Ье 
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<іх г 


ас 
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дЬ° 


ди° 
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( 8 ) 


где к является положительным корнем уравнения 
Ьек 2 +нРк-а=0. 

Решения (6) удовлетворяют граничным услови- 
ям (5) и условиям непрерывности и°(х), Ь°(х) и их 
первых производных в точке х = 0. 

Предел а->со соответствует бесконечно большой 
энергии активации. Тогда и,— >1, и если принять 

1 ~ ц « ~ І 2Т 0 КТ Ь 

и, ~ * Т Ь -Т 0 ^Т Ь Е ’ 


то получим выражение для скорости пламени ѵ„, 
приведенное в [8]. За масштаб скорости ѵ. удобно 
взять ѵ„. Тогда 


и» 0 = 1, а = Ьея 


я -1 Е 

я 777; 
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С учетом сказанного в дальнейшем будем пола- 
гать иА— 1 . 


С принятым видом температуры воспламене- 
ния и. в пределе Е/КТ ь -^<х>, скорость химической 
реакции в форме (2) стремиться, как и закон Арре- 
ниуса, к 5-функции Дирака. 


2. Математическая формулировка задачи для сла- 
бо искривленного пламени. Наложим теперь на 
плоское стационарное пламя малое возмущение, 
вызывающее такую же по величине деформацию 
его фронта по поперечным координатам у, і. Для 
линейного по степеням возмущений анализа дос- 
таточно рассматривать систему уравнений 


^ + хѵ^- = Аи + ІѴ,^ + хѵ^ = ІлАЬ+Ж, 
ді дх ді дх 

Дг зі+ зі+ аі. 

дх 2 ду 2 & 2 ’ 


( 7 ) 


Вычислив входя щ ие в (8) производные 

Щ- = а(1-ъ°)8(и° -и,), ^ = -аг\{и а - и,), 
ди до 

далее имеем 

^ - (О + Х 2 )Р + а{\ -Ь°)8(и° -и.)Р = 

сіх ах 

= ат)(и° -и,)С + <уЯ 2 

Ье ^ ^ - (О + Я 2 Ье)С - сщ(и° - и,)С = ( 9 ) 

= -а( 1 - Ь°)8(и° -и,)Р + а к 2 

сіх 

Граничные условия (5) справедливы и для воз- 
мущенного пламени, только в (5) необходимо про- 
извести замену и°-»и, 5°— >6. Из полученных таким 
образом условий для системы (9) следует обраще- 
ние в ноль Р и С при х — х — го и их первых производ- 
ных при х— >+со. В точке х=<^ и в рамках используе- 
мого здесь линейного анализа должны быть непре- 
рывными значения и, Ь, сіи/ сіх, сіЬ/сіх. Условие неп- 
рерывности и и Ь дает равенства 

р = р, С, = С 2 . (10) 

Условие непрерывности их производных (точ- 
нее, потоков энергии и вещества) приводит к тре- 
бованиям: 

а г и\ + ар, _с ? 2 м° + АР г а 2 Ь/ + ас , _ й 2 Ъ\ + ас 2 ^ л ) 
сіх 2 сіх сіх 2 сіх ’ Дх 1 ах ах 1 ах 

Вторые производные от м°(х), Ь°(х) имеют раз- 
рыв в точке х=0: 

ау__сРиі_ к а 2 ъ ° а 2 ъ\ _ ^ 
ах 2 ах 1 ’ ах 2 ах 2 те ’ 


где пренебрегается конвективными переносами по 
направлениям у, і, как малыми величинами высших 
порядков. Безразмерное время ( измеряется в едини- 
цах к/ и]. Решения уравнений (7) будем искать в виде 

и =и°(х) + ^(у, 2 ,()Р(х), Ъ =Ъ°{х) + ^{у, 2 ,с)С{х), 
хѵ= ѵі'° + ѵі'' = ѵі'° -цА’^, и'' = -^гА'§, 

где в укороченном операторе А' отсутствует диффе- 
ренцирование пох. Деформация Е,{у,і,1 ) фронта горе- 
ния имеет экспоненциальную зависимость от време- 
ни с инкрементом О нарастания возмущений (или 
показателем роста возмущений) и синусоидальную 
зависимость от пространственных координат с вол- 
новыми числами А ь Я 2 по направлениям у и і. 


и тогда условия (11) в окончательном виде прини- 
мают форму 


^±-^± + іс = о с1С ' с, ° 2 I к _ 0 

ах ах ’ ах ах іх 


( 12 ) 


Интегрируя (9) по исчезающей малой области 
вблизи х = 0 и используя известные [10] свойства 
5-функции, находим 


ар 2 

ах 


ар 

ах 


+ {к + 1)77,= О, 


ас, 

ах 


_^_ + к±1 р ■ 0 (1 3) 

ах Ье 


Учитывая (10), легко видеть, что (12) и (13) сов- 
местно приводят к равенствам 

(14) 
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Это означает, что одно из равенств (14) можно 
взять в качестве дополнительного к ( 10 ) и ( 12 ) гра- 
ничного условия. Таким образом, при х=0 имеем 
пять граничных условий на нахождение неизвест- 
ных четырех постоянных интегрирования (после 
выполнении условий при х->±со) и собственных 
чисел ч и О. 

3. Характеристическое уравнение. Приступим те- 
перь к решению уравнений (9). Для области х<0 из 
(9) с учетом ( 6 ) получаем систему уравнений 

і ^- ! *- (П+А!№= '' лг ттт“ р<Л 

и ^^- <п+А ' и,0 '- ,А '^ ехр( й ) - 

Ее решения, обращающиеся в ноль при л:->-оо, 
следующие: 

Рх = Л ехр(ах) ~ Ч я ; к \ і ехр(х), 

1 + л/і + 4(0 + А 2 ) 
а= 2 ’ 

С | .*,ехр(^)-, й ^ ^/ехрО+и), ([5) 

1 + Л /і + 4Ье(0 + ЬеЯ 2 ) 

^ “ 2Іе ' 


Аналогично в области х>0 соответственно имеем 

^ - ( й + 12 ) р 2 = аС 2 + </А 2 ехр(-кх), 

^е^-^-^ 2 --(С^ + ^еX 2 )С^ -аСі = аХ 1 — ехр(-Ьс), 
<Ъс<к 2 2 4 а 


С 2 =ё 2 ехр(-хх)-с/ 


Х = 


П + ЬеЯ- а 

-1 + ^/і + 4Ье(а + П + Я 2 Ье) 
2Іе 


ехр(-іг), 


А =/ 2 ехр(-ух) + А 1§2 ехр(-^х)-^ Л Д, ехр(-Ьт), 

і2 + ЬеЯ 

1 + л / 1 + 4( а + Я 2 ) а (16 

7 2 * 2 +*- ( П + Я 2 )’ 


Ж = 


& Я(Я + 1)-(П + ЬеЯ 2 ) 
к + 1 /Е(Я + 1)-(П + Я 2 ) 


В (15) и (16)/,/, й, & -постоянные интегриро- 
вания. 


Подстановка (15) и (16) в граничные условия 
( 10 ) приводит к алгебраическим уравнениям 




Я 2 к 
П + А 2 к + 1 


/ 2 +г 24 ^а4& Л ’ й=& - (17 > 


а/-? 


Использование (15), (16) в (12) дают уравнения 

к_^„4- 4 х 2 и 2 +к = 0, 

к 


Рёх-д 


п+я 2 к + 1 

Я 2 А : 


+г/тиѴ« — - 

+ Х8і~Ч 


Г2 + ІхА~ 

я 2 я- 3 


П + ЬеЯ 2 аіх 


2 +Г^= (18) 

П + ІхА 2 а Ье 




я 

Ье ’ 


где при записи второго выражения использовано вто- 
рое равенство из (17) и известное выше Ее^ 2 +^-й=0. 
Из второго выражения в (18) находим 


8 1 = 


Я 2 


Ье(я + Р) ѵ О + ЬеЯ 


--1 


(19) 


Из условий (14) определяем/,/: 


/.=■ 


Я 2 


А + Цп + Я 2 


- 1 ,/ 2 = ? 


Я 2 


- 4 &-— ■ ( 20 ) 

П+ЬеЯ 2 А-+1 ' 


Из первых в (17), (18) и (19), (20) после неслож- 
ного расчета получим уравнение 

Я 2 


к 

к + 1 


ЧУ 


-а -у 


Я 2 

О + ЬеЯ 2 


0і = 
02 — 


к Х-У 


о+я 2 

А: 

О + ЬеЯ 
А(ЬеА + 1) 


~02<? ж , г „,2 + * = 0 » 


+ / ж 2 + ж-(0 + Я 2 )’ 
А(А-у) А(А + 1)-(0 + ЬеЯ 2 ) 
к + 1 А(А + 1)-(0 + Я 2 ) 


( 21 ) 


Из условия выполнения требуемого равенства в 
(21) необходимо найти явный вид рО. Параметр 
ц определяет изменение скорости распространения 
пламени при наличии возмущения на первона- 
чально плоском фронте. Для нахождения нового 
значения скорости пламени необходимо разре- 
шить уравнение ( 21 ) при всевозможных значениях 
рО, считая д таким же корнем ( 2 1 ) , как и О . При 
этом в качестве решения задачи ( 1 ) получаем набор 
пар < 7 , 0. Поясним на частном примере, как выпол- 
няется процедура нахождения данных собственных 
значений. Пусть для простоты Ье=1. Тогда 


л, = Х~У к\к + 1) 

% Х+Р Г + ж-(П + Я 2 ) Х + Р’ 


02 


к(к-у) 
к + 1 


и вместо ( 21 ) имеем 


--1 


О + Я 2 


Х-У 
Х + Р 


■ТЗ 2 '-* 


Равенство в (22) выполняется при 
« = Л + ,, 


( 22 ) 

(23) 


и таких значений О, которые являются корнями 
уравнения 


Х-У к ~ 2 У = п 
Х + Р к + і 


Решая это уравнение, находим 0,—Я 2 , 
0,=-Я 2 -1 /4. Подставив их в (23), получим #(О|)=0, 
<7(0 2 )=-1/(4Я 2 ). 

С учетом формулы (23) определение скорости 
пламени у>=мР-дА'І; теперь принимает форму 

— °-(іі+ 1 ] л ^. 


Заменим здесь О и Я 2 на их выражения через 
производные от деформации фронта согласно ра- 
венствам 


П 1 ,2 1 д,Е 0 д,Е 

П = ?ЗГ’ Я = -^:у>=у> +-ЙГ-А5- 


Первый член #=1 - это скорость плоского не- 
возмущенного пламени, второй - кинематическая 
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добавка, появление которой обусловлено тем, что 
задача решается относительно системы координат, 
связанной с плоским невозмущенным пламенем, 
третий определяет изменение скорости пламени от 
кривизны его фронта. 

Вопрос о диффузионно-тепловой устойчивости 
ламинарного пламени можно сформулировать дву- 
мя разными способами. В первом из них скорость 
ѵ„ распространения пламени можно считать не ме- 
няющейся при наличии возмущений. Тогда при 
< 7=0 устойчивость определяется отсутствием в кор- 
нях (21) инкремента О с положительной действи- 
тельной частью. Такой способ применен в работах 
[4, 5,7]. Второй же подход, использованный выше, 
опирается на предположении непостоянства ско- 
рости пламени (<р^0) в возмущенном пламени. Ка- 
кой из них адекватно описывает реальное пламя, 
предстоит еще выяснить. Но ясно, что в первом 
способе накладывается ограничение на спектр воз- 
можных значений О. 

4. Устойчивость горения при постоянной скорости 
пламени. Анализ устойчивости пламени при 
ѵ„=сошІ(<7=0) изучен достаточно подробно [8]. По- 
этому данный вопрос здесь разбирается только по- 
верхностно. Считая число Льюиса произвольным, 
рассмотрим (21) в пределе к->оо. Полученное та- 
ким образом уравнение 

</=Г°-+Ье1 Ѵ^-Р + У 

и 2 Л/Ье-/3+у+А 2 (Ье-1)/а 

при <7=0 распадается на два других: Г2/А 2 +Ье=0; 
1/Ье -/1+у=0, откуда находим 0,=-ЬеА 2 , 0 2 =Я. 
Второй, положительный, корень указывает на аб- 
солютную неустойчивость пламени в случае беско- 
нечно большой энергии активации. Этот результат 
получен ранее в работе [7], но следует также и из 
дисперсионного соотношения [4] . 

При <7=0 из (21) получим 

к -2у х~У Аг(ЬеАг + 1) _ р 

к + 1 Х + Р Ье(Ѵ+х-(0 + А 2 )] 

Следует заметить, что это характеристическое 
уравнение при к» I не переходит в аналогичные 
выражения, полученные в [4, 7]. 

В работе [5] было выявлено, что при числе Ье>1 
существуют области с колебательной и неколеба- 
тельной потерей устойчивости. Граница области 
для первого случая находится из условия 0=0, для 
второго - из условия 0=/ш, /=-1. Параметр со име- 
ет смысл частоты. Если же Ье<1, то наблюдается 
только колебательная потеря устойчивости. Но ис- 
следование устойчивости, проведенное на основе 
формулы (25), показало несколько другую картину. 
Так, при Ье>1 наблюдается только неколебатель- 
ная (апериодическая) потеря устойчивости. Ре- 
зультаты численного анализа (25) приведены на 
рисунке, а. Здесь области неустойчивого горения, 
занимающие пространство выше кривых 1 , 2 и т.д. 
при различных значениях безразмерной энергии 
активации к имеют форму полуостровов с общим 


берегом - ординатой Ье. Если же Ье< I , то потеря 
устойчивости носит только колебательный харак- 
тер. Причем область неустойчивого горения имеет 
форму острова, размеры которого быстро растут с 
увеличением к. Результаты численного расчета 
приведены на рисунке, б. 




Рисунок. Изменение области апериодической (а) потери ус- 
тойчивости (1.е>1) с ростом безразмерной энергии 
активации к: 5 (1), 10 (2), 15 (3), 20 (4), 30 (5); пери- 
одической (б) потери устойчивости (іе<1) с ростом 
безразмерной энергии активации к: 5 (1), 7 (2). При 
к=5 остров неустойчивости образован кривой 1 и 
частью оси іе, а при к=7, 10, 15 еще и частью оси А 2 

5. Устойчивость горения при переменной скорости 
пламени. Найденные выше для случая Ье= 1 показа- 
тели Оі=-А 2 , ГУ=-А 2 -1/4 роста возмущений с от- 
рицательно определенными знаками, казалось бы, 
свидетельствуют об абсолютной устойчивости пла- 
мени с переменной ѵ„. Тем не менее, оказывается, 
что кроме ГУ и ГУ существует еще один режим го- 
рения с инкрементом ГУ, который выпадает из 
внимания, если в (21) сразу полагать Ье=1. 
Действительно, пусть в уравнении (21 ) 
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Естественные науки 


я 2 


■ Ье. 


В получающемся из (21) равенстве 
. О + ЬеЯ 2 


У 


-а - у -(к -у)х 


П + Я 2 

х М± 1)-(^ + ^ 2 ) + я +1 = о 

Я(Я + 1)-(0 + Я 2 ) 


(26) 


примем 0.=С1)=-Х 2 +к(к+ 1). Раскрывая возникаю- 
щую при этом неопределенность вида 0/0 в дроби 
к -у 

к(к + 1)-(0 3 +Я 2 )’ 
приведем (26) к виду 

=0. (27) 

(к + 1)(2к + 1) ѵ ’ 

Требуемое здесь равенство нулю левой части 
(27) возможно только при Ье=1, если не принимать 


во внимание экзотический случай Я=0, соответ- 
ствующий одномерным возмущениям и нереаль- 
ный случай к= со. 

При инкременте 0 3 в решениях Т, и Р 2 тоже воз- 
никает неопределенность вида 0/0, которая легко 
раскрывается. Громоздкий, но простой расчет дает 
случай Ье= 1 : 

х<0, Р, = С, = --г^— ехрх; 

1 1 к + 1 

х > О, Р 2 =С 2 = - + ехр(-Ах). 

Таким образом, для волновых чисел Я, удовлет- 
воряющих неравенству 

Л<ущт 1)»*=^, 
пламя неустойчиво. 
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